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The performance of telescope systems working at microwave or visible/IR 
wavelengths is typically described in terms of different parameters according 
to the wavelength range. Most commercial ray tracing packages have been 
specifically designed for use with visible/IR systems and thus, though very 
flexible and sophisticated, do not provide the appropriate parameters to fully 
describe microwave antennas, and thus to compare with specifications. In this 
work we demonstrate that the Strehl ratio is equal to the phase efficiency 
when the apodization factor is taken into account. The phase efficiency 
is the most critical contribution to the aperture efficiency of an antenna, 
and the most difficult parameter to optimize during the telescope design. 
The equivalence between the Strehl ratio and the phase efficiency gives the 
designer /user of the telescope the opportunity to use the faster commercial 
ray-tracing software to optimize the design. We also discuss the results of 
several tests performed to check the validity of this relationship that we 
carried out using a ray-tracing software, ZEMAX and a full Physical Optics 
software, GRASP9.3, applied to three different telescope designs that span a 
factor of ~ 10 in terms oi D/X. The maximum measured discrepancy between 
phase efficiency and Strehl ratio varies between ~ 0.4% and 1.9% up to an 
offset angle of > 40 beams, depending on the optical configuration, but it is 
always less than 0.5% where the Strehl ratio is > 0.95. 
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1. Introduction 

Performance evaluation is a critical step in the design of any optical system, either at 
microwave or visible/IR wavelengths. The image quality criteria more commonly used, 
however, are quite different in these two regions of the electromagnetic spectrum. In 
fact, in the analysis of microwave antennas and radio telescopes the two fundamen- 
tals figures-of-merit used by designers and users are the aperture efficiency and the 
beam efficiency, whereas in optical systems the Strehl ratio and ray aberrations are 
often quoted. This is because of the coherent nature of most microwave antennas, 
where single-moded receivers are generally used (exceptions may be millimeter and 
submillimeter bolometers used in radio astronomy), making the phase distribution in 
the image as important as the amplitude distribution in determining the performance 
of the optics. In fact, the aperture efficiency is intrinsically dependent on the phase 
distributions since it is calculated as a correlation integral between the focal region 
field produced by an incident plane wave and the horn aperture field. 

The difference between the microwave and visible/IR wavelengths regimes, in terms 
of the image quality criteria applied to astronomical telescopes, has been reduced 
over the past 10-15 years thanks to the development of focal plane arrays (FPA, 
hereafter). In fact, the noise performance of receivers used in radio astronomy has 
improved dramatically during this time, especially at millimeter and submillimeter 
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wavelengths. As a consequence, it has become clear that the best means of increasing 
observing efficiency for mapping extended sources or to conduct blind searches is to 
use imaging arrays located at the focal plane of the telescope. This implies the need 
of a larger field of view (FOV) with few aberrations in the range of frequencies used 
by the array(s) of receivers. Very often these FPA require some relay optics to convert 
the telescope focal ratio (which, in some cases, may be quite large, i.e. ^ 10) to the 
smaller focal ratios of the individual feed-horns. As a consequence, the overall image 
quality of the total system, telescope and reimaging optics, must be evaluated over a 
wide FOV, thus effectively contributing to bridging the gap between the microwave 
and visible wavelengths regimes. 

A number of commercial ray tracing packages exist that are being used to analyse 
the performance of FPAs for use with existing or planned (sub) millimeter telescopes. 
However, many of these packages have been specifically designed for use with optical 
(i.e., visible and IR) systems and thus, although very flexible and sophisticated, they 
do not provide the appropriate parameters to fully describe microwave antennas, and 
thus to compare with specifications. The possibility to easily convert an optical-based 
design parameter, such as the Strehl ratio, to a fundamental antenna-based design 
parameter, such as the phase efficiency, gives the designer/user of the telescope the 
opportunity to use the faster commercial ray-tracing software to optimize the design. 
Once the design is optimized, a full Physical Optics software can be used to analyse 
more thoroughly all critical performance parameters of the antenna (e.g., spillover, 
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antenna noise temperature, etc.). Another advantage offered by this conversion con- 
sists of the possibihty to study tlie degrading effects on the wavefront caused by 
obstructions to the beam (e.g., secondary reffector and its support struts) which are 
notoriously difficult to simulate in Physical Optics software. 

In this paper we review the main design parameters generally used in evaluating 
the performance of optical designs at both microwave and visible wavelengths. Based 
on this review we find a simple relationship between the (antenna-based) aperture 
efficiency and the Strehl ratio. We also show the results of several tests performed to 
check the validity of this relationship that we carried out using a ray-tracing software, 
ZEMAX and a full Physical Optics software, GRASP9.3, applied to three different 
telescope designs. 

The paper is organized as follows: In Sect. [2] we review and discuss the definitions 
of antenna gain and aperture efficiency while in Sect. Owe analyse the definition of 
Strehl ratio and derive a simple relationship between the aperture efficiency and the 
Strehl ratio; in Sect. IHwe show the results of a comparison obtained using a Physical 
Optics and a ray-tracing program and, finally, we draw our conclusions in Sect. O 

2. Antenna gain and aperture efficiency 

2. A. Definitions 

The gain of an antenna is a measure of the coupling of the antenna to a plane wave 
field, and it can be written in terms of the effective area (we assume that ohmic losses 
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are negligible): 



For an aperture type antenna the gain is expressible in terms of the illumination 
by the feed. We can assume that the illumination is linearly polarized, and that the 
aperture lies on an infinite plane. In this case the gain is expressible in terms of Ea{r'), 
the magnitude of the (in-phase) illuminating electric field in the aperture plane. If 
almost all of the energy in the field is contained in a small angular region about the 
z' axis, and if we use the scalar-field approximation, then G{6, 0) can be written as:^ 
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4:71 



^a(r',R) dS' 



AP' 



JE^^{v')dS' 



(2) 



with 



£a(r',R) = Ea{v')e^''^'">e^^^-'' (3) 
R ■ r' = r' sin 6 cos(0 — 0') 
dS' = r'dr'd(t)' 

where we have introduced the complex electric field in the aperture, i'a(r',R). We 
have also indicated with k = 27i/X the wavenumber, and the field point Q at position 
r' on the aperture plane (see Fig.[T]) has polar coordinates (r', 0'). R is the unit vector 
along the direction to the observation point, with 6 representing the angle formed by 
the direction to the observation point and the optical axis and being the angle 
measured in the plane of scan, i.e. perpendicular to the optical axis (z'), as shown 



in Fig. [H The integral in the numerator is calculated over the antenna aperture, 
whereas the integral in the denominator must extend over the entire plane if there is 
any spillover illumination in the case of reflector antennas. 




Fig. 1. Coordinate systems used to calculate the antenna gain. 



The phase aberration function, $(r'), in Eq. defines the phase at point r' in the 
aperture plane, which accounts for any change in the optical path length resulting 
from the structural deformation of the primary reflector, the displacements of the 
secondary reflector and the feed. Thus, it is in $(r') that one can take into account 
the positions of different feed-horns in a FPA. 
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For aperture type antennas, the effective aperture can be related directly to the 
antenna geometric area, Ag, by means of the aperture efficiency, rj^{6,(f)) (e.g., see 



Ref. 2) 



A, 



(4) 



Therefore, 



G{9, 



^a(r',R) dS' 



AP' 



AJEl{v')dS' 



(5) 



(6) 



The on-axis gain. Go, is obtained by setting R ■ r' = 0, then we obtain: 



G„ 



4:7lA, 



A2 



AP' 



(7) 



AgJE^Ar')dS' 



If the phase is constant over the aperture the on-axis gain attains its maximum value, 



G 



G 



M 



^2 



Ea{r') dS' 



AP' 



AJEl{v')dS' 



(9) 
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A case of special interest is that of uniform illumination over the aperture, i.e. 



Ea{r') = const over the antenna aperture and zero outside. Hence, we obtain 77^^ = 1 



and the ideal gain, Gidoah is then defined as 



AnA, 



(11) 



M • 



Thus, we obtain the well-known result that the uniform field distribution over the 
aperture gives the highest gain of all constant-phase distributions over the aperture.^ 

2.B. Phase- error effects 

In the previous section we showed that if the phase distribution is constant over the 
aperture, the maximum gain, Gu, is obtained in the direction of the optical axis, i.e. 
R- r' =0. However, if a phase-error distribution is present over the aperture, this may 
no longer be the case. A phase-error over the aperture, i.e. deviations from uniform 
phase, may arise from various causes, such as a displacement of the feed-horn from 
the on-axis focus (e.g., in FPAs), or distortion of the optical surfaces, or it may be 
caused by phase-error in the field of the feed-horn. 

If the phase distribution is a linear function of the aperture coordinates, then it 
can be shown that the far-field is the same as that of the constant-phase distribution 
but displaced with respect to the z'— axis, i.e. the direction of peak-gain is no longer 
in the direction of the system optical axis.^ In the case of arbitrary phase ditributions 
over the aperture, if the phase-error does not deviate too widely from constant phase 
over the aperture, and if it can be decomposed into a linear phase distribution and 
higher-order terms, then we may write 




(12) 
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Fig. 2. Off-axis feed and tilted aperture plane geometry. 
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where $i(r') is linear in the coordinates over the aperture and causes an undistorted 
beam shift, i.e. a change in direction of the peak gain (now corresponding to 6' = ^pk), 
whereas $ab(r') accounts for the true wave front distortion. The shifted far-field beam 
can then be considered to have arisen from a tilted aperture plane, i.e., from the 
aperture projected onto a plane normal to the direction of the peak gain, Rpk, as 
shown in Fig. O In the projected aperture the linear phase distribution term cancels 
out, leaving only higher-order phase errors, i.e.. 



^a(r, R) dS 



AP 



(13) 



A,JE^,{T)dS ' 

oo 

where now 

£:a(r,R) = E,(r)e^■*-(^)e^■'=^■^ (14) 

and where r is the position of a point in the projected aperture plane, indicated with 
AP, such that R ■ r = for R = Rpk- If ^'ab(r) = then the field distribution has 
constant phase over the projected aperture and the antenna gain in this aperture will 
be given by,^ 

Gup = Gm cos6'pk (15) 

where cos 6'pk — 1 for most radio astronomical applications. Therefore, in the following 
sections we will refer to the antenna gain and aperture efficiency as the gain and 
aperture efficiency in the projected aperture plane, unless noted otherwise. 



11 



2.C. Main contributions to the aperture efficiency 

The aperture efficiency of an antenna is determined by a number of phenomena and 
hence it can be written as the product of a number of individual contributions (e.g., 
see Ref. 3): 

r]^{9, (p) = ?7spill ?7taper(6', 0) VphaseiO, 0) (16) 

where ?7spiii is the spillover efficiency, ?7taper is the taper efficiency and ?7phase takes 
into account all phase-error effects causing a distortion of the wave front. We have 
also assumed that ohmic losses are negligible and that the aperture is unblocked. 
The spillover efficiency includes all spillover contributions from the feed, subreflector, 
diffiaction, etc., 



El{r)dS 

Vsm=-r • (17) 



oo 



?7tapcr accounts for the aperture illumination taper due to the feed and the reflector 
geometry. 



'7taper ( 



AP 



2 



A,fEl{v)dS ' 

AP 



and finally, ?7phasc accounts for the residual high-order phase distortions of the wave- 
front at the aperture plane, due to optical aberrations, surface errors or misalignments. 
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etc., 



^a(r, R) dS 



AP 



AP 



2 • 



(19) 



In the direction of the peak gain R ■ r = 0, as we earher mentioned, and thus the 
R = {9, (p) dependence can be dropped from r/taper and ?7phase- 

In the case of on-axis, dual-reflector systems the central subreflector and its support 
structure cause a partial shadowing of the aperture, which leads to a loss of efficiency. 
To take this effect into account the integral at the numerator of Eq. ( |T3l) can be 
written in the case of a partially blocked aperture: 



£^{v,K)dS= l£^{v,K)dS- I S,{r,R)dS (20) 



APu 



AP 



subr 



where APuock represents the area of the aperture plane subtracted of the blocked 
part, AP indicates as usual the full area of the aperture plane and subr indicates the 
integration area over the subreflector, assuming this is the main source of blockage. 



By substituting Eq. ( 1201) into Eq. ( 1131) we thus obtain 



^a(r, R) dS 



AP 



A,!El{v)dS 



X 



^a(r,R) dS 



^a(r, R) dS 



AP 



(21) 



where the first term at the right can once again be written as in Eq. (fT6l) and thus 
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the second term can be interpreted as the blocking efficiency due to the subreflector, 



AP 

We note that in the direction of the peak-gain (Rpk ■ r = 0), for an uniform, unaber- 
rated ($ab(r) = 0) field we find the well-known result, 



where Aprim and Asuhr are the surface areas of the primary and secondary reflectors, 
respectively. In general, the geometrical blockage caused by the support struts can be 
up to several times larger than the blockage caused by the secondary mirror, especially 
in open-air antennas. Therefore, the blockage efficiency given by Eq. ( l23ll usually 
overestimates the real efficiency and should be corrected including the strip blockage 
of the plane-wave and the blockage from the converging spherical-wave between the 
primary mirror and the subreflector (e.g., see Ref. 4). 

3. Strehl ratio 

3. A. Strehl ratio on- axis 

While the main antenna-based figures-of-merit are usually, though not necessarily, 
defined in the far field of the aperture, the Strehl ratio of an optical imaging system is 
defined as the ratio of the aberrated to unaberrated incoherent Point Spread Function 
(PSF, hereafter^). 




(22) 




(23) 
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Fig. 3. Coordinate frame at exit pupil {x',y', z') and position, p', of point P 
at focal plane. The field amplitude at point r' on the system's exit pupil is 
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When considering the optical system in receiving mode, the PSF refers to the in- 
stantaneous field distribution in the focal plane of an optical imaging system produced 
by a far-field point source. For simplicity we assume that the fields are emerging from 
the exit pupil of the optical system with a system focal length /, and converging 
towards the image plane. Let's suppose that the exit pupil is on an infinite plane 
located at 2; = 0, and with the normal unit vector in direction of the z-axis, n = z' 
(see Fig. [3]|j. Then, following Ref. 5,6, the scalar field at a point P at position p' in 
the paraxial focal plane (see Fig. [3]) is given by 

Efip') oc J E,,{v')e-^7^'''' dS' (24) 

AP' 

where E^.^ is the field amplitude at a point Q at position r' on the system's exit pupil 
and / is also equal to the radius of curvature of the reference sphere centered at point 
O in the focal plane. In the case of a point source in the far field of the system Eex 
is uniform over the pupil. 

In Eq. (I2II) the substitution of the exit pupil for the antenna aperture plane, and 
the consequent use of r' in both cases, is justified by using the equivalent parabola 
(e.g., in a dual-refiector system) and by the fact that when the point source object is 
at infinity, then the diameter (assuming a circular aperture) of the exit pupil can be 
substituted with the diameter of the entrance pupil, or main dish in a dual-refiector 



^The focal plane and the observation point in the far-field defined by R in the previous sections 
(where the optical system was considered in transmission mode) lay on opposite directions with 
respect to the x'y' plane. This will be taken into account in Sect. 13. CI 

16 



system (see Ref. 6, p. 184), and the system focal length would be in this case the focal 
length of the equivalent parabola/ In other words, the spherical (i.e., aberration-free) 
wavefront leaving the equivalent parabola and converging to the focus is identified 
here with the Gaussian reference sphere centred on the exit pupil. Then, we can state 
that the (unaherrated) incoherent PSF is simply the square modulus of Ef{p'), i.e., 



PSF = /(pO = \Efip 



/M2 



Eq. fl24j) is strictly valid in the absence of phase errors that may modify the per- 



fectly spherical convergent wave that was assumed earlier in the special case of an 



aberration-free wave- front. In the more general case of a distorted wave- front Eq. (l24l 



should be re-written as: 



Ef{p') oc I E,,{v')e-'l'''-'' e^^'^^'US' 

AP' 



(25) 



where $(r') is the phase error term. The Strehl ratio, S, of the imaging system is 
then given by the ratio of the central (i.e., p' = 0) irradiance of its aberrated and 
unaberrated PSFs. From Eq. ( 125|) S can be written in the form^'^: 



So 



m 



J Ee.(r')e^*(^')d5' 

AP' 



m\ 



$=0 



/ E,,{r')dS' 



AP' 



(26) 



where 5*0 = S{p' = 0). The Strehl ratio can also be used as a measure of the on-axis 
PSF away from its central irradiance peak, and thus we can write: 



S{P') 



I{P') 



$=0 
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(r') e-^7'''-^'e^*('-') dS' 
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J E,,{r')dS' 








AP' 







(27) 



3.B. Strehl ratio off-axis 

In equations (I24p to (1271) the position in the paraxial focal plane of the central irra- 
diance peak of the PSF was taken as the origin of a Cartesian system of axes and 
also as the center of the (unaberrated) Gaussian reference sphere.® The observation 
of an object point off-axis, which is equivalent to having the feed lateraly displaced in 
a microwave antenna, introduces both a change in the position of the PSF peak (or 
direction of peak gain in an antenna) and wave-front aberration. The quasi-spherical 
(i.e., aberrated) wave will be thus converging to a point displaced with respect to 
point O in Fig. [31 If p'p^ represents the position of the off-axis PSF peak in the focal 
plane, then Eq. (1271) can be re-written as: 



S{P') 



Hp') 

Hp'pk)\'S>,b=0 



AP' 



J Eex(r')e-^7'''pk"'e^*i(''')rf^' 



AP' 
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where <I>ab and $i(r') have been defined in Eq. f|T2|) and Sect. 12. B] Thus, /(p'pk)|$ab=o 
represents the peak irradiance of the nnaberrated, off-axis PSF. 

In Sect. I2.BI we saw that by tilting the aperture plane so that it becomes perpen- 



dicular to the direction of the peak gain, it is possible to write the aperture efficiency 



in terms of $ab only. Likewise, in the definition of the PSF it is possible to align the 
2;-axis along the direction from the center of the exit pupil to the off-axis Gaussian 
image point, which can also be taken as the origin of a new Cartesian system of axes. 
The Gaussian image point is also the center of curvature of the (tilted) wave front, 
and for this point all path lengths from the spherical wave front would be equal, in 



the absence of higher-order phase distortions. Then, Eq. ( l28l) takes the same form as 



Eq. (I27D, i.e 



Sip) 



AP 



AP 



(29) 



where the peak of the PSF is now at point p = in the new system of axes, centered 
on the Gaussian image point in the focal plane, and r now lies on a tilted plane, AP, 
perpendicular to the direction of the off-axis PSF peak. Thus we have in the projected 



plane. 



So 



/(O) 



/(o)l 



*ab=0 



/ ^ea:(r)e^*-Wt/5 
AP 



AP 



(30) 



3. C. Strehl ratio and aperture efficiency 



In this section we use the previous results to derive a relationship between aperture 



efficiency and Strehl ratio. First, we use Eq. (1131) to form the ratio of the aberrated 
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and unaberrated aperture efficiency (in the projected aperture plane), i.e. 



with 



where ?7^(R) = f]^ 



direction R 



S.{r, R) dS 



AP 



Ea{v) dS 



AP 



(31) 



Ea{r) dS 



AP 



(32) 



AJEl{v)dS 

oo 

0) is the aberrated aperture efficiency measured in the generic 
[i.e., not coincident with the direction of the peak gain, Rpk), 



for the general case in which the direction of peak-gain is not along the main optical 
axis of the system, as explained in Sect. 12. B[ rj^^ = r7^(Rpk)|$^,^=o is the unaberrated 
aperture efficiency measured in the direction of the (ofT-axis) peak gain, i.e. r/MP 
represents the peak aperture efficiency as measured in the projected aperture plane. 
Recalling that in the direction of the peak-gain Rpk ■ r = (see Sect. 12. CI) the In- 
dependence can be dropped from ?7mp. From eqs. ( ITOl) and ( |T5l) it also follows that. 



(33) 



if << 1, where 77,^ is the maximum aperture efficiency as defined in Sect. 12. A[ 
From Eq. (!32|) and equations (fT7|) and (fT8|) we also see that rj^ = r^spiu r/taper- 

Then, we note that equations fl29p and (1311) have the same form and, for small 



angles close to the optical axis it holds that 



a ■ r = — R • r 
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where we have defined a = p// (see the discussion in Ref. 9). However, since E^x 
represents the field produced by a point source in the far field of the system, in order 
to conclude that equations fl29l) and fl3Tl) are fully equivalent one must assume that 
the incident field on the optical system from a distant source has an apodization 
equivalent to that produced by the feed illumination on the antenna aperture (see 
Sect. KM . In this case we can write -E'ex(i') = -E'a(r), and thus 



(34) 



Then, by comparing equations (fT6l) to (fT9l) with Eq. (l34ll one can see that in general. 



7]^ S{p) = ?7spill?7taper(R.)?7phase(R-) 



(35) 



Usually, however, one is interested in the aperture efficiency at the nominal position 
of the peak gain (i.e., at the center of the far- field beam), or equivalently at the center 



of the PSF, then it also holds that 



— ''7m ^° 



Vm ^spill ^taper 



(36) 



and 



'S'o — ''^pha 



(37) 



with Rpk • r = and r]o = VaI^ ~ R-pk) is the aperture efficiency in the direction of 
the peak-gain, corresponding to Eq. ([8]) in the projected aperture plane, i.e. 



J Fair) e^^-i^W dS 

AP 



A,SEl{v)dS 

oo 
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(3^ 



where we have not used the subscript "p" (for "projected parameter") in rjo because 
of the approximation in Eq. (l33ll . Therefore, Eq. [3l\ ) finally shows the equivalence 
between the Strehl ratio and phase efficiency. 

Clearly, r/^ takes into account both taper and spillover effects, whereas is a 
measure of the phase aberrations. Therefore, in the case of an nnaberrated wave 
front, i.e. 5*0 = ?7phasc = 1, the aperture efficiency is rjo = r]^ and depends only on the 
spatial distribution of the field over the antenna aperture. Furthermore, by explicitly 
writing the aberration function, <l>ab(i")j in terms of the primary aberrations (e.g., see 
Ref. 5) it would be possible to derive the individual contributions to the aperture 
efficiency by, e.g., coma, astigmatism and curvature of field, which are usually the 
most relevant aberrations in radiotelescopes. However, this is beyond the scopes of 
this work and will not be done here. 

4. Comparison of Strehl ratio and aperture efficiency 

In this section we want to compare the values of the Strehl ratio, obtained from a 
ray-tracing optical software, ZEMAX (Focus Software^"), and the associated value of 
?7phase5 obtained through the numerical integration of Eq. (fT9|) and using the aperture 
field values computed by a Physical Optics program, GRASP9.3 (TICRA Engineering 
Consultants^^). Several configurations have been analysed and are discussed below. 
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4- A. Description of software packages 

The analysis has been conducted using the GRASP9.3 package, which is a com- 
mercial tool for calculating the electromagnetic radiation from systems consisting 
of multiple reflectors with several feeds and feed arrays. This package can use sev- 
eral high-frequency techniques for the analysis of large reflector antennas, such as 
Physical Optics (PO) supplemented with the Physical Theory of Diffraction (PTD), 
Geometrical Optics (GO) and Uniform Geometrical Theory of Diffraction (GTD), 
which require a moderate computational effort. 

The PO technique is an accurate method that gives an approximation to the sur- 
face currents valid for perfectly conducting scatterers which are large in terms of 
wavelengths. The PO approximation assumes that the current in a specific point on 
a curved but perfectly conducting scatterer is the same as the current on an infinite 
planar surface, tangent to the scattering surface. For a curved surface, the PO current 
is a good approximation to the actual one if the dimensions of the scattering surface 
and its radius of curvature are sufficiently large measured in wavelengths. The well- 
known GO method uses ray-tracing techniques for describing wave propagation. Since 
GO gives discontinuities in the total electromagnetic field, GTD is often applied in 
addition to GO, since GTD methods may account for diffraction effects. 

On the other hand, ZEMAX is a classical optical design tool based on ray-tracing 
methods, which combines three major categories of analysis in one package: lens 
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design, physical optics, and non-sequential illumination/ stray light analysis. 

4-3. Calculation of the aperture efficiency with GRASP9.3 

As described in Sect. 14. At GRASP9.3 allows several methods for the electromagnetic 
analysis of the reflecting surfaces. An interesting tool of GRASP9.3, based on the 
ray-tracing, for calculating the aperture field is the so-called "Surface Grid".^^ This 
method returns the reflected magnetic field on the surface according to the formula: 
Hr = Hi — 2n(n ■ Hi), where Hi is the magnetic incident field and n is the normal 
to the surface. The magnetic reflected field on the surface, Hj., is then projected, 
with a phase adjustment, on the aperture plane. As described in Sect. 12. Bl when 
the feed is placed off-axis the aperture plane is tilted according to the direction of 
the peak-gain. For a dual reflector configuration, the scattering from the secondary 
and primary mirrors has been analyzed through the GTD technique and the "Surface 
Grid", respectively. 

This approach is particularly appropriate when the diameter-to-wavelength ratio 
of the primary reflector is very large and when the observation point is in the near- 
field (such as the aperture plane case). Under these conditions the PO method would 
be very time-consuming; in fact, it would require a huge number of points on the 
reflector where currents need to be evaluated. Using the method described here to 
analyze the primary reflector the diffracted field from the edge of the reflector is 
not considered. However, the numerical results obtained with this "hybrid" technique 
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have been compared with those obtained by applying the PO method to both the 
primary and secondary mirrors, resulting in a very good agreement between the two 
methods. 

In order to calculate the aperture efficiency from Eq. ( 13T|) we use the the complex 
electric field in the aperture plane, i.e. £^a(r, R), produced by GRASP9.3, which is 
tabulated through its real and imaginary components. These can then be used to 
calculate the amplitude and the phase function of the field. The complex electric field 
is finally read by a proprietary code which evaluates Eq. ( JTOl) in order to determine 
the phase efficiency. 

4-0. Comparison of results 

The values of the Strehl ratio and phase efficiency obtained with ZEMAX and 
GRASP9.3, respectively, have been compared using three different optical systems. 
These systems have been selected to represent standard telescope designs, and the 
frequencies used in the simulations cover the mm- and submm-wavelength regimes. 
For the electromagnetic analysis with GRASP9.3, we have always used a linearly po- 
larized Gaussian feed. Although more realistic feed models to describe circular horns 
could be adopted, for the sake of comparison with ZEMAX and to avoid introducing 
any systematic error due to different feed illumination, we report the results obtained 
with a Gaussian model only. The level of apodization in ZEMAX has then been chosen 
to be consistent with that produced by the Gaussian feed-horn. 
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4.C.I. Single-dish antenna 

First, we have carried out the comparison in the simplest possible case, i.e. an un- 
blocked spherical reflector antenna. This choice eliminates or minimises potential 
discrepancies due to different handling in ZEMAX and GRASP9.3 of effects such as 
multiple reflections, aperture blocking and diffraction at secondary surfaces. 

95,0 -1 , , , , , , , 1 




88,0 > 

87,0 J ^ ^ ^ ^ ^ ^ ^ , 

0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 

Beam tilt [deg] 

Fig. 4. Plot of the Strehl ratio and of the phase efficiency at a wavelength of 
500/im for the case of a spherical reflector 105 cm in diameter with a //# = 2. 

The surface chosen for this simulation is spherical because it ensures that spherical 
aberration will limit the overall FOV to small 1°) angles near the optical axis. This 
is required in order to avoid introducing further variables in the comparison between 
ZEMAX and GRASP9.3 due to the incidence angle of radiation over the aperture of 
the feed-horn in the focal plane, which may affect the coupling between the PSF and 
the electric fields on the horn aperture. The selected aperture was 105 cm in diameter 
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with a. f = 2 and the simulations have been carried out at a wavelength of 500 
fiBi. For the electromagnetic analysis with GRASP9.3, a linearly polarized Gaussian 
feed has been used with a taper level of —12 dB at 14°. 

The results are shown in Fig. HI the comparison has been extended up to a maximum 
offset angle of ~ 1.4°, or about 44 beams at 500 /xm, and the maximum measured 
difference between the Strehl ratio calculated by ZEMAX and the phase efficiency 
calculated by GRASP9.3 is 0.38% at the maximum offset angle. We also note, however, 
a 0.25% discrepancy on boresight, which will be discussed in the next section. 

4.C.2. Dual-reflector antenna: Cassegrain configuration 

We have then analysed the most common radio telescope design, consisting of a dual- 
reflector antenna. We flrst consider the classical Cassegrain conflguration, which we 
have derived from the design of the "Balloon-borne Large Aperture Submillimeter 
Telescope" (BLAST) telescope. Compared to the original design with a spherical 
primary mirror^^ and to the newer telescope design with a Ritchey- Chretien optical 
conflguration, the system analysed here has a parabolic primary and a hyperbolic 
secondary. The diameters of primary and secondary mirrors are 181.61 and 42.76 cm, 
respectively, and the system focal ratio is 5. As in the single- reflector linearly 
polarized Gaussian feed has been used, but with a taper level of —9 dB at 6°. 

The results are shown in Fig. (5) the comparison has been extended up to a maxi- 
mum offset angle of ^ 0.79°, or about 42 beams at 500 /im, thus quite equivalent to 
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Fig. 5. Plot of the Strehl ratio and of the phase efficiency at a wavelength 
of 500/im for the case of a classical Cassegrain telescope. The diameters of 
primary and secondary mirrors are 181.61 and 42.76 cm, respectively, and the 
system focal ratio is 5. 
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Fig. 6. Same as Fig. [5] for the scaled-up version of the BLAST telescope. The 
primary and secondary reflector diameters equal to 12.2m and 2.6m, respec- 
tively. 
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the previous simulation. The maximum measured difference between the Strehl ratio 
calculated by ZEMAX and the phase efficiency calculated by GRASP9.3 is about 
0.59% at an offset angle of about 0.5°. We observe that the discrepancy between the 
two methods is also relevant (0.2 — 0.3%) for offset angles near boresight and it is 
possibly more systematic in this case than in the single-reflector design analysed in 
the previous section. 

This on-axis difference is likely due to the relatively small secondary diameter to 
wavelength ratio, -Dgcc/A, which may cause an on-axis decrease of the antenna gain 
due to diffraction effects from the edge of the secondary. To test this hypothesis, 
we have scaled-up the BLAST telescope, while keeping constant the wavelength, in 
order to obtain an optical design with a much larger -Dsec/ A ratio, comparable to that 
used in the next section for the "Sardinia Radio Telescope". We have thus obtained 
a telescope with the same focal ratio at the Cassegrain focus but with a primary and 
secondary reflector diameter equal to 12.2m and 2.6m, respectively. 

The results are shown in Fig. (6) in this case the comparison has been extended up 
to a maximum offset angle of ~ 0.33°, or about 116 beams at 500 /im. As expected, 
the discrepancy near the optical axis has decreased compared to both the single-dish 
and the original BLAST cases. The maximum difference is about 0.61%, thus still 
quite similar to that observed in the original BLAST design despite the much larger 
offset angle in beam units used in the scaled-up telescope. These results indicate that 
diffraction effects are calculated differently in GRASP9.3 and ZEMAX. 
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4.C.3. Dual-reflector antenna: Gregorian configuration 

The third system analysed during this comparison is another dual-reflector antenna, 
though in a Gregorian conflguration. In this case we have changed the wavelength 
to a larger value of 3 mm and have also chosen a telescope with a much higher D /\ 
ratio. The baseline design is in this case the "Sardinia Radio Telescope" (SRT^^); 
however, we have converted the original shaped design of the SRT to a more standard 
Gregorian conflguration, keeping the same aperture (64 m) and system focal ratio 
(2.34) of the SRT. As in the previous two cases, a linearly polarized Gaussian feed 
has been used, with a taper level of —12 dB at 12°. 

The results are shown in Fig. [T) the comparison has been extended up to a maximum 
offset angle of ^ 0.136°, or about 42 beams at A = 3 mm, thus consistent with the 
simulations used for the single-dish and the BLAST conflgurations. The maximum 
measured difference between the Strehl ratio calculated by ZEMAX and the phase 
efficiency calculated by GRASP9.3 is about 1.9%, thus larger than in the optical 
systems discussed above. However, in the range of offset angles where the Strehl 
ratio (or equivalently the phase efficiency) is > 0.95, i.e. the range which is normally 
targeted by the optical design of diffraction-limited telescopes, the difference between 
Strehl ratio and phase efficiency is < 0.5%, consistent with that observed in the 
BLAST telescope. 
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Fig. 7. Plot of the Strehl ratio and of the phase efficiency at a wavelength of 
3 mm for the case of a classical Gregorian telescope, with a primary reflector 
diameter of 64 m and a system focal ratio of 2.34. 

5. Conclusions 

We have reviewed the main design parameters generally used in evaluating the per- 
formance of optical designs at both microwave and visible wavelengths. In particular, 
we have reviewed the classical concept of antenna gain and the main contributions 
to the aperture efficiency, with special attention to phase-error effects. We have then 
described the formalism with which to compare the aperture efficiency and its com- 
ponents with the Strehl ratio, which is the standard parameter used to evaluate the 
image quality of diffraction-limited telescopes at visible/IR wavelengths. We have 
shown that a simple relationship can be found between Strehl ratio and aperture ef- 
ficiency: the Strehl ratio is equal to the phase efficiency when the apodization factor 
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is taken into account. 

We have then compared these two parameters by running ray-tracing software, 
ZEMAX and full Physical Optics software, GRASP9.3, on three different telescope 
designs: a single spherical reflector, a Cassegrain telescope and finally a Gregorian 
telescope. These three configurations span a factor of ^ 10 in terms of D/\. The 
simple spherical reflector allows the most direct comparison between Strehl ratio and 
phase efficiency, as it is only marginally affected by edge diffraction effects. In this 
case we find that these two parameters differ by less than 0.4% in our ZEMAX and 
GRASP9.3 simulations, up to an angle of about 44 beams off-axis. The other two 
configurations are more prone to diffraction effects caused by the secondary reflector, 
especially in the case of the smaller Cassegrain telescope. 

The phase-efficiency is the most critical contribution to the aperture efficiency of the 
antenna, and the most difficult parameter to optimize during the telescope design. The 
equivalence between the Strehl ratio and the phase efficiency gives the designer /user 
of the telescope the opportunity to use the faster (and less expensive) commercial 
ray-tracing software to optimize the design using their built-in optimization routines. 
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